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The response of vortex state to the magnetic field in NbsSn is probed using muon spin rotation 
and small-angle neutron scattering. A transformation of vortex structure between hexagonal and 
squared lattice is observed over a relatively low field range of 2-3 Tesla. The gradual increase of the 
magnetic penetration depth with increasing field provides microscopic evidence for anisotropic (or 
multi-gapped) s-wave pairing suggested by the Raman scattering experiment. This result renders 
need for careful examination on the difference of electronic properties between NbaSn and V3SL 

PACS numbers: 74.25.Jb, 74.25.Nf, 76.75.+i, 61.12.Ex 



While the technologies for utilizing high-T c cuprates 
are making steady progress in recent years, triniobium 
stannide (NbaSn) is still one of the most important ma- 
terials for the practical application of superconductivity 
after fifty years since its discovery.! Because of its conve- 
nient characters such as modestly high transition temper- 
ature (T c ~ 18.3 K) and remarkably high upper critical 
field (fioH C 2 ~ 24.5 T), NbaSn is used at the heart of 
superconducting devices, most notably as superconduct- 
ing wires in high- field magnets. NbaSn belongs to a class 
of A%B binary intermetallic compounds with the A15 or 
/^-tungsten structure, where those based on Nb give rise 
to a subclass that includes Nl^Ge^, which used to domi- 
nate the highest T c (~ 23 K) for more than thirty years, 
and therefore it was once a subject of intensive study 
until '70s AAA Surprisingly, however, it seems that the 
microscopic details on the superconducting characters of 
NbaSn are not fully understood to date. For example, the 
revelation of a structural transformation from cubic to 
tetragonal with decreasing temperature not far above T c 
(martensitic phase transition, which occurs at Tm = 43 
K in NbsSn£) led to the argument that the coupling to 
the instability of lattice vibrations (soft phonons) associ- 
ated with the martensitic transition might contribute to 
the increase of T c . Experimental verification of the elas- 
tic softening 7,8 prompted various models for microscopic 
mechanism of the martensitic transition, many of which 
emphasize the importance of the quasi-one-dimensional 
character of Nb atoms in the A15 structure that might 
also lead to the anomaly in the electronic density of 
states, the Kohn anomaly due to electron-phonon inter- 
action, or to the Peierls instability^ 4 *^ Unfortunately, 
while these models were partially successful to explain 
the elastic properties of NbaSn, none of them has ap- 
peared as unambiguous explanation for the high T c . 

Here, we report on the first systematic study of vor- 
tex (flux line) state in NbaSn using muon spin rotation 
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FIG. 1: a) Example of /xSR time spectrum (noH — 2 T) 
displayed in the rotating reference frame of 201 MHz (where 
solid curve represents fitting result using a model described 
in the text, yielding a reduced chi-square of 1.21), and b) real 
amplitude of the fast Fourier transform (FFT, raw data) of the 
time spectra under several magnetic fields in the mixed state 
of NbaSn (oc n(B)), where each spectrum is shifted vertically 
for clarity. The background signal from muons stopping in a 
lucite scintillator (used both as a sample holder and a muon 
veto counter) gives rise to the peak at B — fi^H — 0, where 
the amplitude is determined by the efficiency of the muon 
veto counter (~ 95 %, which leads to a relative background 
yield of -5%). 



(/iSR) and small-angle neutron scattering (SANS). The 
combined use of /iSR and SANS on the identical spec- 
imen is proven to be ideal for the study of the mixed 
state, as they provide complementary information on the 
magnetic field distribution associated with flux line lat- 
tice (FLL): /iSR is sensitive to local field distribution 
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while SANS is sensitive to the long range structure of 
FLL. Moreover, it eliminates much of the problems asso- 
ciated with the difference in the preparation of specimen. 
In this brief communication, we show that one of the 
keys to understand the superconductivity of NbaSn is the 
anisotropy in the electronic state, which is demonstrated 
by the transformation of FLL structure from hexagonal 
to squared with increasing flux density. Similar observa- 
tion in VsSi strongly suggests that this is a common fea- 
ture of the superconductors with A15 structure jliiSL where 
the anisotropy of the Fermi surface is reflected through 
the nonlocal effect Meanwhile, the field-induced en- 
hancement of the magnetic penetration depth at lower 
fields, which is absent in VsSi, 13 indicates the presence of 
low energy quasiparticle excitation due to the quasiclassi- 
cal Doppler shift. 1 — This result strongly suggests that the 
superconducting order parameter also has an anisotropic 
(or multi-gapped) structure on the Fermi surface™ 

NbaSn samples were prepared by the method described 
elsewhere. 16 The superconducting transition tempera- 
ture was determined to be 17.9(2) K, as inferred from the 
shielding effect in magnetization measurement. Several 
single crystals with a sizable dimension (~ 10 2 mg each, 
covering a beam spot area of 7x 7 mm 2 with random crys- 
tal orientation) were used for /iSR measurement with a 
high transverse field spectrometer (HiTime) installed on 
the M15 beamline of TRIUMF, Canada, where they were 
irradiated by a muon beam of 4 MeV. Further details on 
the /iSR measurements using HiTime might be found in 
earlier publications™ One of these crystals, which was 
confirmed to be of single domain at ambient temperature, 
was used for the neutron measurements with the SANS- 
J spectrometer of JRR-3, Japan Atomic Energy Agency. 
In both measurements, precautions were taken to mini- 
mize the effect of flux pinning: every time when external 
magnetic field was changed, the sample temperature was 
raised above T c and cooled down again after the magnetic 
field was settled. The magnetic field dependence was ob- 
tained for transverse field /iSR spectrum (TF-/iSR, where 
the initial muon spin polarization is perpendicular to the 
external field) at 2.0-2.5 K and for SANS spectrum at 
2.6(1) K, respectively. 

Since we can reasonably assume that implanted 
muons (with 4 MeV) stop randomly on the length 
scale of FLL (which is 20-200 nm over the relevant 
field range), TF-/iSR signal provides a random sam- 
pling of the internal field distribution B(r), P(t) = 
exp(— (Jpt 2 ) n(B) exp(i7 /i 5t — i(j))dB„ with n{B) = 
(S(B — B(r))) r , where <j p is the additional relaxation due 
to random flux pinning, n(B) is the spectral density for 
the internal field defined as a spatial average (( ) r ) of the 
delta function over a unit cell of FLL, 7 M is the muon 
gyromagnetic ratio (= 2tt x 135.53 MHz/T), and <j) is 
the initial phase of rotation™ These equations indicate 
that the real amplitude of the Fourier transformed muon 
spin precession signal corresponds to the spectral den- 
sity, n(B). In the modified London model, which is a 
good approximation to evaluate B(r) at relatively lower 
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FIG. 2: SANS diffraction patterns obtained by subtracting 
the background data at 30 K from the data taken at 2.6(1) 
K after field cooling for an external magnetic field fioH = 1 
T (a), 2 T (b), 3 T (c), and 4 T (d) (raw data). The crys- 
tal c-axis is aligned to neutron beam, and one of the {1,0,0} 
in-plane axis is aligned to the horizontal axis in the plane 
perpendicular to the field. 



external fields (fioH <C i^ C 2) 19 , B(r) is given by a sum of 
the magnetic induction from isolated vortices to yield 
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where K are the vortex reciprocal lattice vectors, Bo 
(~ IJLqH) is the average internal field, A is the effective 
London penetration depth, ^ v is the cutoff parameter for 
the magnetic field distribution near the vortex cores, and 
rixxyy and d are dimensionless parameters arising from 
the nonlocal corrections with the term proportional to 
K^Ky controlling the fourfold anisotropy: the overline 
denotes an averaging of the parameter over the polar an- 
gle of the symmetry axis, as it is anticipated for a mosaic 
of crystals. 11 It is evident in Eq.QJ that b(K) is a spa- 
tial Fourier component of B{r). The model predicts an 
asymmetric field profile for n(B) characterized by a neg- 
atively shifted sharp cusp due to the van Hove singularity 
associated with the saddle points of B(r) and an asymp- 
totic tail towards higher fields where the maximal field is 
determined by B(\r — ~ £ v ) (with being the vortex 
centers). 

Examples of the fast Fourier transform (FFT) of /iSR 
spectra at several external magnetic fields are shown in 
Fig. d where one can readily observe the asymmetric 
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lineshape with slight smearing due to the random pin- 
ning of vortices, random local fields from nuclear mo- 
ments, and the limited time window for FFT (~ 4 /is: 
note that such smearing is irrelevant for the analysis in 
time domain). Based on the least-square method with 
appropriate consideration of the statistical uncertainty, 
the fiSR spectra are compared in time domain with those 
calculated by using Eqs.QJ-Q to deduce a set of param- 
eters, <j p , A, fy, n xxyy , and d: a typical example is shown 
in Fig. H For the reconstruction of B(r) from n(B), 
one needs to know the two-dimensional structure of FLL 
obtained from other experimental techniques including 
SANS. In particular, it is noticeable in Fig. ^ that the 
peak due to the saddle points exhibits a steep change 
with increasing external field from 2 T to 3 T, suggesting 
an anomaly in the FLL structure. As shown below, this 
speculation is confirmed by SANS measurements over the 
region of relevant magnetic field. 

In a SANS experiment, a FLL gives rise to Bragg re- 
flections at reciprocal lattice points, K = lu+mv = K^ m . 
The intensity of a single reflection (Z,m) is given by 
|6(K)| 2 /|K/ m |. Considering Eq.©, it is reduced to 7 10 oc 
dio^o/A 4 f° r a fi rs t order reflection, where K\$ = 2it / d\Q 
with dio being the FLL interval (= y/3ao/2 for a hexag- 
onal FLL) and $o is the flux quantum. Fig. [2| shows the 
typical SANS- J data at 2.6(1) K after subtracting back- 
ground data obtained at 30 K (> T c , which is far below 
Tm), where scattered neutrons from an incident beam 
with a mean wavelength Ao = 6.5 A (AAo/Ao —0.13 
FWHM, collimated within a diameter of 2 mm) are 
detected by a position-sensitive detector (consisting of 
128x128 pixels with each pixel size being 4.95 x 5.08 
mm 2 , covering a sensitive area of 58 cm ) at a distance of 
10 m (1-2 T) or 4 m (3-4 T) from the specimen. One of 
the crystal c-axes and external magnetic field are aligned 
together to the neutron beam. The cryostat (including 
the magnet) and specimen are rotated slightly (±0.2°, in 
several steps) around the vertical axis to obtain a sum of 
diffraction patterns over Bragg angles for both triangu- 
lar and square FLL structure. It is clear in Fig. that 
there is a difference in the pattern of angles between 
data at 1-2 T and those at 3-4 T. Since the patterns 
at lower fields suggests an overlap of two (possibly dis- 
torted) hexagonal domains due to the tetragonal crystal 
structure, their change into the nearly four-fold symmet- 
ric pattern at higher fields demonstrate the occurrence 
of triangular-to-square lattice transformation at 2-3 T. 
The field range of transformation is considerably higher 
than that in VaSi, where the transition to square lat- 
tice is nearly complete at 1.5 T (with H ||{0,0,1})£ As 
shown in Fig. 0Jl, the apex angle of the real space unit 
cell is evaluated to be indeed close to 60° and 90° for the 
respective data. Furthermore, the spots on the {1,0,0} 
axis at high fields (Figs. d) is absent at lower fields. 
This clearly indicates that the diagonal of the real space 
rhombic unit cell, which is aligned to {1,1,0} or {1,1,0} 
direction in each domain, rotates to a {1,0,0} direction 
with increasing field. A similar domain structure at lower 



field has been observed in VsSiA2£ 
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FIG. 3: Spectral density distribution n(B) calculated for a 
hexagonal FLL with (d — 0.004) or without (d = 0) non-local 
correction under an external field of 0.5 T (where A = 83.5 
nm, f„ = 3.1 nm, and n xxyy = 0). The sharp peaks are due 
to van Hove singularities arising from turning points in B(r), 
where V5(r) = 0. An additional peak (labeled by arrow) ap- 
pears for the former case because of the incommensurability 
between triangular FLL and squared shape of B(r) around a 
vortex. Those for a square FLL at 4 T are shown for compar- 
ison (A = 100.0 nm, £ v — 2.8 nm, and n xxyy — 0). 

Once the FLL structure is resolved by SANS measure- 
ments, high quality /iSR spectra can provide rich infor- 
mation on the mixed state though the detailed map- 
ping of B(r); as it has been demonstrated in the ear- 
lier works 13, — , the lineshape strongly depends on A, 
n XX yy and d in Eq.Q and thereby they can be deduced 
by fitting analysis without much ambiguity. In particu- 
lar, as shown in Fig. |3 the spectral weight for the lower 
field side of the van Hove singularity strongly depends 
on the magnitude of d for the square FLL. Besides these, 
we need to introduce two more parameters to describe 
the reciprocal vectors which undergo transformation: the 
apex angle, /?, and the mean width for the FLL distor- 
tion, <tx, 22 to describe the fluctuation of the reciprocal 
vectors, K = Z(l + 5j)u + ra(l + 5 m )v with a probability 
distribution function, P(Sk) oc exp(— ^/cr 2 K ) (k = i,m). 
The last parameter effectively serves to add a Lorentzian- 
type broadening to n(B) and thereby discernible from the 
Gaussian broadening due to local disorder (which is de- 
scribed by cr p ). The randomness associated with the mo- 
saic crystals is presumed to be represented by ctk- Fig.0] 
shows the field dependence of those parameters, where 
(a) A, (b) fy, (e) n xxyy and d are directly related with the 
local structure of supercurrent around a single flux line, 
whereas others are with the morphology of the FLL. It is 
inferred from the behavior of gk (Fig- St) that a strong 
long-range distortion occurs over a field range of 2-4 T, 
where the FLL structure undergoes transformation be- 
tween triangular and square lattice. The anomalous en- 
hancement of A observed in Fig.^i is probably due to the 
artifact associated with this distortion; the reduced chi- 
square exceeds 1.5 around 3 T while it is typically 1.2-1.4 
for fji H < 2 T or fi H >4 T (see Fig.gf). Similar scat- 
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tering of the parameter value is also noticeable in £ v , a p , 
and d over the relevant field range, which may represent 
the systematic uncertainty due to the distortion. How- 
ever, except for these anomalies, we can observe a clear 
trend in all of the parameters. Namely, A exhibits a grad- 
ual increase with increasing magnetic field below ~ 2 T 
(which we discuss below in more detail), while ^ v tends to 
decrease slightly at lower fields and levels off above ~ 2 T. 
The behavior of £ v , which to some extent represents that 
of the effective vortex core radius, is commonly found 
in various superconductors, which is attributed to the 
stronger overlapping of supercurrents around vortex cen- 
ter and intervortex transfer of quasiparticles. 23 The mag- 
nitude of at lower field is in good agreement with the 
Ginzburg-Landau coherence length, £gl — 39 A deduced 
from the relation H C 2 = ^o/2tt^q L [with H C 2(2 K) being 
assumed to be 22 T]. As is anticipated from the SANS re- 
sult, (3 exhibits a steep change from 60° to 90° around 2-3 
T. The behavior is qualitatively in good agreement with 
the prediction of the non-local model (solid curv o 11 ! 12 , 
where the critical field is assumed to be 3 T). Thus, it 
indicates that the anisotropic flow of the supercurrent 
due to the nonlocal effect is the primary origin of the 
square-lattice formation at higher fields. The magnitude 
of d (~ 4 x 10 -3 ), which should be regarded as a lower 
boundary considering the fact that it is reduced by an av- 
eraging over a varying symmetry axis^i is considerably 
larger than that observed in VsSi (~ 0.3-1 x 10 -3 ), 13 sug- 
gesting a stronger anisotropy in the Fermi surface. The 
spin relaxation due to the random local distortion (a p ) 
exhibits a gradual decrease with increasing field: this is 
qualitatively understood as a result of stronger overlap 
of flux lines at higher field that makes less room for flux 
lines to deviate away from the regular position. Unfortu- 
nately, the SANS data are available only up to 5 T due to 
experimental conditions, and therefore the origin of the 
weak anomaly around 5.5 T (see Fig.^h) is not identified 
at this stage. 

It is Interesting to note that the anisotropy parame- 
ter d is sharply reduced to zero below ~2 T: a similar 
behavior is reported for V3SL 13 As shown in Fig. EF, a 
comparative analysis with d fixed to 0.004 (~ an average 
for hqH > 2 T) yields systematic increase of chi-squares, 
while little variation is found for other parameters. This 
is attributed to an additional cusp in n(B) (see Fig. |3| 
corresponding to the saddle points relatively near the vor- 
tex cores where B(r) is influenced by the squared shape 
of supercurrent flow. Our result indicates that the field 
profile over the field range of hexagonal FLL may be bet- 
ter described by the isotropic London model. 

Now, we discuss the origin of gradual increase of A in 
Fig. Our systematic study on the magnetic field de- 
pendence of A in various superconductors indicates that 
the gradient (dX/dH) provides information on the low en- 
ergy quasiparticle excitations^ In the mixed state, the 
quasiparticle momentum vp is shifted by the flow of su- 
percurrent v s around the vortices due to a semi-classical 
Doppler shift, leading to a shift of the quasiparticle en- 
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FIG. 4: Magnetic field dependence of various parameters to 
describe the field distribution B(r) obtained by analyzing TF- 
/iSR spectra using the modified London model: (a) the effec- 
tive London penetration depth (A), (b) the cutoff parameter 
(£u), (c) the enhanced relaxation rate due to local disorder of 
FLL (<7 P ), and the mean width for the long-range distortion 
((7k), (d) the apex angle of the rhombic unit cell in the real 
space (together with the SANS data), and (e) the dimension- 
less parameters to describe anisotropy of supercurrent density 
around the flux lines (n xxyy , d). Error-bars in (a)-(e) show 
statistical errors (those not shown are smaller than the sym- 
bol size), whereas the scatter of points, particularly for 2-4 T, 
may represent systematic uncertainty (see text). The reduced 
chi-squares are shown in (f). The open symbols in (a)-(d) jmd 
(f) (only for fioH < 2 T) are those obtained by fixing d to 
0.004. Dashed lines are merely intended to guide for eyes. 
[See text for the solid lines in (a) and (d).] 



ergy spectrum to an amount e = ravp • v s . When the su- 
perconducting order parameter has a nodal structure or 
anisotropy (including multi-gap), the Cooper pairs with 
a gap energy of less than e can be broken by the Doppler 
shift, 14,24 leading to an enhancement of A. Using a nor- 
malized magnetic field, h = H/H C 2, we can express the 
field dependence of A, 



X(h) = A(0)[l+Tyfc], 



(3) 



where rj is the dimensionless parameter to represent the 
gradient (the magnitude of pair breaking effect). Since 
X~ 2 (h) = 47re 2 n s (/i)/m*c 2 that is determined by the su- 
perfluid density, n 3 (h), the reduction of n s (h) due to the 
Doppler shift would lead to the increase of X(h). A fit- 
ting analysis of the data in Fig. ^1 by Eq.© below 2 T 
yields rj = 1.24(2) with A(0) = 818(1) A, where the re- 
sult is shown by a solid line in Fig. Considering that 
no clear evidence is reported for the presence of nodes 
in the order parameter, the situation in NbaSn is close 
to the case of MgE>2 where a similar trend of X(h) with 
r] c± 1.3 was observed. 17 Since /iSR measurements are 
made at finite temperature (T), the Cooper pairs hav- 
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ing the energy near the gap minimum (or a smaller gap) 
A s would be susceptible to excitation due to the Doppler 
shift when A s < /c#T. 24 Interestingly, an earlier report 
on electronic Raman scattering experiment suggests a rel- 
atively strong gap anisotropy (~ 20%) in NbaSn, whereas 
no such anomaly seems to exist for VsSi^ Meanwhile, a 
very recent specific heat measurement in NbaSn demon- 
strated the presence of residual excitation below ~ 0.3T C , 
which they attributed to a double-gapped order param- 
eter with A s = 0.61(5) meV™ Our result is another, 
yet microscopic evidence for the anisotropic (or multi- 
gapped) s-wave pairing for the order parameter of NbaSn. 
The anisotropy of the Fermi surface revealed by FLL 
transformation may be related with the structure of the 
order parameter. Further study using the techniques sen- 



sitive to gap anisotropy would be helpful to determine the 
gap structure in full detail. 

In summary, we have obtained detailed mapping of the 
local magnetic field in the mixed state of NbaSn using 
/iSR and SANS, from which we extracted fundamental 
length scales of superconductivity. Unlike VsSi, the field- 
induced enhancement of the effective penetration depth 
strongly suggests that the order parameter in NbaSn has 
a structure, either anisotropic or multi-gapped. 
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